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ABSTRACT 
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O , We consider the effects of rapid pair creation by an intense pulse of 7-rays propagat- 

ing ahead of a relativistic shock. Side-scattered photons colliding with the main 7-ray 
beam amplify the density of scattering charges. The acceleration rate of the pair-loaded 
pL^ ' medium is calculated, and its limiting bulk Lorentz factor related to the spectrum and 

compactness of the photon source. One obtains, as a result, a definite prediction for 
the relative inertia in baryons and pairs. The deceleration of a relativistic shock in the 
{Sj . moving medium, and the resulting synchrotron emissivity, are compared with existing 

results for a static medium. The radiative efficiency is increased dramatically by pair 
loading. When the initial ambient density exceeds a critical value, the scattering depth 
tranversed by the main 7-ray pulse rises above unity, and the pulse is broadened. This 
O ' sets an upper limit to the pre-burst mass loss rate of ~ 1CP 5 M per year, and places 

significant constraints on gamma-ray burst progenitors. An anisotropic 7-ray flux (on 
an angular scale T _1 or larger) drives a large velocity shear that greatly increases the 
energy in the seed magnetic field forward of the propagating shock. 

6 

Subject headings: gamma-rays: bursts - theory - radiation mechanisms 



1. Introduction 

Physical models for 7-ray emission from a relativistic fireball, and the ensuing synchrotron 
emission from the decelerating shock, have generally neglected the feedback of the intense 7-ray 
flux on the dynamics of the fireball, or on the prompt and delayed emission. There are, however, 
several reasons to believe that this interaction can have an important influence on the observed 
multiwavelength emission from gamma-ray bursts (GRBs). Pair creation by high energy photons 
will raise the radiative efficiency of a shock (Thompson 1997); and - if the ambient medium is 
sufficiently dense - it will significantly smear and broaden the observed time profile of the prompt 
GRB. The pairs also carry net forward momentum which will couple effectively to the ambient 
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medium if this contains a strong enough magnetic field. In this work we extend previous results on 
radiative acceleration in the Thomson (Noerdlinger 1974; O'Dell 1981; Phinney 1982; Kovner 1984) 
and Klein-Nishina (Madau & Thompson 1999, hereafter Paper I) scattering regimes to include the 
effects of pair creation. Attenuation of the high energy 7-ray spectrum by collisions with soft 
photons was calculated by Baring and Harding (1997), but the dynamical effects of pair creation 
and the feedback of radiative acceleration on the scattering depth were not considered in that work. 
The role of pair creation in accelerating a trans-relativistic outflow from a black-hole accretion disk 
with a hard X-ray spectrum has been considered recently by Beloborodov (1999). 

The plan of this paper is as follows. In § 2 we calculate the rate of acceleration by a thin shell 
of 7-rays, some of which side-scatter off ambient charges and collide with the main 7-ray beam. 
We show that acceleration is much faster (for a given compactness) than in the Thomson limit, 
because the highest energy photons provide most of the momentum. The reduction in the mean 
inertia per scattering accompanying pair creation significantly increases the limiting Lorentz factor 
over the value for a baryonic plasma. In § 3 we apply these results to GRB afterglows. We consider 
the deceleration of a relativistic shock propagating into a medium that is itself moving at bulk 
relativistic speed, and show how the time scaling of the synchrotron emission is modified from the 
case of a static medium. We estimate the limiting ambient density above which the 7-ray pulse 
experiences a large optical depth to scattering. Finally, we note that angular variations in the 
radiative force will excite non-radial shear flows in the ambient medium that can greatly amplify 
the pre-existing seed magnetic field. 

In the following we will denote by x and x s the incident and scattered photon energies in units 
of m e c 2 in the (unprimed) lab frame, and with T = (1 — /3 2 ) -1 / 2 the bulk Lorentz factor of the flow. 

2. Pair loading 

Gamma-ray sources of high compactness are opaque to photons that move at an angle with 
respect to the 7-ray flow. When the source spectrum extends to an energy x>l, even soft photons 
will create pairs via 77 — > e + e~ , with a cross-section close to Thomson near threshold. This implies 
that the source of 7-ray radiation must itself expand at relativistic speed. When the photons are 
beamed into a narrow angle along the direction of motion, the threshold energy for pair production 
within the beam is increased to x ~ In Paper I we calculated the radiative acceleration of 

matter exposed to a strong source of 7-rays, including the full Klein-Nishina suppression of the 
scattering cross-section, but did not include the effects of pair creation. However, the same Compton 
scattering process that accelerates material sitting close to a 7-ray source also creates side-moving 
photons. When the integrated flux of photons of energy x > 1 is large, each side-scattered photon 
deposits its entire momentum, along with the momentum of a second colliding photon, into an 
electron-positron pair. The build-up of a pair plasma ('pair loading') through this process, and the 
bulk acceleration that ensues, will be the subject of this section. 
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2.1. Radiative Acceleration 

At very high compactness, the source of gamma-rays must be impulsive; otherwise the as- 
sumption of low scattering optical depth in the accelerating medium is not self-consistent. We 
approximate the photon source as a radially propagating shell situated at c(t — At) < r < ct. 
We found in Paper I that, in the Thomson limit, a test particle surfing the photon shell attains a 
limiting Lorentz factor 

r max = 5(^ (i) 



2 V A* 

when the shell is thin enough that the acceleration length is much less than the radius, 

frLxcA* « r. (2) 

The parameter 

r c = o^At — — , (3) 
is expressed in terms of the integrated energy flux 

(4) 



J F x dx, 



and characterizes the column density of photons iV 7 through the shell, r c ~ otN 1 {x). When the 
particles are coupled by a magnetic field and the electrons and positrons are cold, the mean mass 
per scattering charge is 

fi = m e + m p ^1 + 2^^ . (5) 

(We assume for simplicity a hydrogen plasma and neglect the inertia of the original neutralizing 
electron component, but not of the pairs.) The term m e in this expression must be replaced with 
the relativistic inertia fi e = ^('y e )m e when the electrons and positrons are hot. 

The optical depth to pair creation is proportional to the optical depth parameter r c , which 
is in turn proportional to the duration At of the gamma-ray pulse, and the source compactness 
Lot I (47r \m? R) . Here, R and L are the source radius and luminosity, respectively. We assume that 
the incident spectrum is broadband, extending well above and below x ~ 1. As long as r c S> 1, even 
a soft photon that is side-scattered can be absorbed by the main photon beam. The accelerating 
material is, in this situation, photon rich: a side-scattered photon that finds itself above the pair 
creation threshold has a much higher probability of colliding with another photon than of Compton 
scattering a second time. Indeed, the number of scattering particles in the relativistically expanding 
material that emits the photons is suppressed by a factor of with respect to the number of 
photons. 

We now calculate the momentum deposited when a (soft) photon B is side-scattered through 
an angle 9 S and collides with a (hard) photon A. The combined process of scattering followed 
by pair creation deposits a net forward momentum (x A + x B )m e c in the accelerating medium. 
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When the scattering charge is cold, the initial energy of the scattered photon in the rest frame 
of the scattering charge is x' = xT(l — (3). A similar relation, x' s = x s T{l — f3cos9' s ), holds for 
the photon after scattering. The scattering angle in the frame of the charge is given by cos9' s = 
(cos 9 S — — (3 cos S ). The lab frame energy x B before scattering can be related to the energy 

xf following scattering through 

B = Txf (1 - /3cos 9 S ) 
X {l-(3)T-(l-cos9 s )xf , {> 

which approaches x B — > xf(l — j3 cos 9 S )/{1 — (3) in the Thomson limit (x' = x' s ). In these variables, 
the cross-section for the collision between photons A and B (scattered) is 



3 



(3 - £ 4 ) In ( 1±|) - 2£(2 - e 



(7) 



(Jauch & Rohrlich 1976), where 



y x A xf (1 — cos 9 S ) ' 
The cross-section reaches a maximum of 0.25(7^ at x A = 2x A lin (^ = 0.71). 
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Most of the momentum deposited through this process involves the side-scattering of very soft 
photons, which collide with hard gamma-rays to produce energetic (and almost radially moving) 
pairs. The minimum (threshold) energy of the radial photon can be expressed as a function of the 
energy and propagation angle of the side-scattered photon, 

A B ^ 
Xmin (Xs ' ° s } = xf(l-cos^) ■ (9) 

Thus, a photon of energy xf <^ (x) that is side-scattered through an angle ~ 9 S creates a pair of 
energy x pair ~ 4:/9^x B . During the first stages of the acceleration, this pair will be relativistic in 
the bulk frame, 

^.pair 

2 7e (injected) ~ - — » 1. (10) 

(Recall that the photon beam is assumed to be radially streaming.) 

In the presence of a sufficiently strong magnetic field, this momentum is effectively commu- 
nicated to the other charges. Inspection of eq. (37) in Paper I shows that this coupling rapidly 
becomes stronger as T increases. In the case of a 7-ray fireball of duration At ~ 10s and photon 
optical depth t c ~ 300 (eq. ||), propagating into an interstellar medium with Bq ~ 3 x 10~ 6 G, 
the coupling is strong only for pairs of lab frame energy x pmr ~ few. However, above a Lorentz 
factor 

Bo V 1/s /T c /At\ 1/5 (xV**\ 2/S 



3 ' 5 U^G 35^ W <»> 
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even energetic pairs deposit their entire momentum. In the process, all particles will share a bulk 
radial Lorentz factor. In the calculations that follow, we assume that -Bo is strong enough (and At 
long enough) to assure effective coupling; but this constraint should be kept in mind when applying 
the formulae below. 

When t c 3> 1, each newly created e + e~ Compton cools on the main photon beam. Indeed, 
it drops to an energy far below the injection energy over the time that the medium accelerates; 
this starting assumption will be justified a posteriori in §2.2. In the process of scattering off an 



anisotropic radiation field, the cooling particle absorbs an additional momentum - similarly to what 
happens in the Compton rocket (O'Dell 1981; Phinney 1982). The effect of this Compton afterburn 
on the net momentum deposited in the medium depends on whether the Compton-scattered photons 
are absorbed by pair creation. When they are not absorbed, the net effect is to increate the total 
momentum absorbed per photon collision to 



/ x (x A + x B )m e c. 



(12) 



Because the magnetic field becomes nearly transverse as the medium accelerates, the numerical 
factor works out to / = | (Paper I). Secondary pair creation can be neglected after the bulk 
Lorentz factor reaches a critical value, but during the early stages of the acceleration the side- 
scattered photons are absorbed. The net effect then is to cancel the additional momentum absorbed 
by the cooling charge, and to force / — > 1. We discuss further details of the pair-photon cascade 



induced by the injection of an energetic pair in §2.3, and for now summarize this process in terms 
of mean relativistic inertia fj, e 



^('j e )m e of the daughter e , and a multiplicity factor iV mu i t 
secondary pair creation can be neglected, one has -/V mu i t ~ 1; otherwise 



2<7e> Wmult 



x A + x B 



When 



(13) 



l + (3)V 

The quantity on the RHS is the total energy of the injected pair in the bulk frame. 

Let us now calculate the rate of acceleration due to injection of energetic pairs which cool 
down rapidly off the main photon beam, to a mean Lorentz factor (jg) -C (x A + x B ). In the fluid 
approximation the rate of change of the bulk momentum is given by 



Of 



r (3c(2fj, e n' e+ + m p n' ) 



+ V 



f3T 2 Pc 2 {2n e n' e+ +m p n' p ) 



J(x A + x B )f 



F x {x J 



x- 



,.B 



-<r 77 (l — cos6* s ) dx A dx B dVt. 



(14) 



where the bulk-frame charge densities are denoted by a prime. The rate pair creation induced by 
side-scattering of soft photons per unit volume, 



d_ 
di 



(T< + ) + V • (Tf3cn' e+ ) = f N ma / x{x ^ )I{x *^ 



-cr 77 (l — cos9 s ) dx A dxf dQ. 



(15) 



The Compton force acting on transrelativistic or subrelativistic particles has been neglected. 
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Substituting the equations of continuity into ( |l4| ) yields 
dT/3 dF(3 



m e (2fi e n e + + m p n p )c 5 



[ \{x A + x B )f - N mnlt 2Tp\ I{Xa B a) °n (l-cos0 s ) dx A dx B dfi. 

(16) 



The value of the second term in the brackets depends on the importance of secondary pair creation. 
During the last stages of acceleration (/ ~ |) one has fi e ~ m e and iV mu it ~ 1 (§[^)- Otherwise, we 
substitute / = 1 and expression fli~3|) for iV mu it , and find that the quantity in brackets is {x A + x B ) 
multiplied by (1 — + f3). Both these results can be expressed in terms of a /3-dependent 

Compton afterburn parameter 



5 



(no secondary pair creation); 



^ \ 1 \+p (secondary pair creation). ^ ^ 

The photon shell will be parameterized by a radial coordinate w, measured from its inner 
boundary, 

r = c(t- At) + w. (18) 

The scattering occurs at position w s and time t s , and the collision at position w± and time t±. 
Hence, 

w s — w± = c(t± — t s ) cos S . (19) 

The intensity l(xf,0 s ) of side-scattered photons is related to the flux F x (x B ) of unscattered radi- 
ation in the thin shell (cf. eq. 0) via 



rcAt fa r , 

I(xf ,0 a ,w±) = J (2n e + + n p ) F x (x B ) (1 - /#)— exp -r 77 (xf , S , w s - m±) dw s . (20) 



Note that the radial position w s of the scattering site within the photon shell always lies outside 
the position w± of pair creation, because the radial speed of the side-scattered photon is less than 
c. The optical depth to pair creation experienced by the side-scattered photon is 

t~Jx b ,0 s ,zu s -w±) = ( [ dx A ——j - (1 ~ cos0 s )a 1 Jx A /x A in (xf ,0 S )] . (21) 

Jvj s C Jx A . (xBfi s ) X A L J 

When the optical depth along the ray from w± to w s is much larger than unity, the momentum 
deposited by side-scattering and pair creation is localized, and the x^-integral can be performed 
directly: 

j t (T(3) = M> J {{x A ) + x B ) F x (x B )(l - P) da^§-, (22) 
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where we have made use of eq. (|5|) for the mean mass fi per scattering charge. The average energy 
of the undeflected, radial photon is 

/ A\ = I F x (x A )a 1 ~ t dx A 

after weighting by cross-section and photon-spectrum. We write 

(x A }(xf,e s ) = k(a)x A Lhl (xf,0), (24) 

from here on. The constant k is plotted as a function of high energy spectra index in Fig. 1. 

It is now necessary to prescribe the radiation spectrum. We assume a simple broken power-law 
form, appropriate for GRB sources, with photon number flux (per logarithm of energy) constant 
below a break energy Xb r ~ 1 and photon index a above the break: 

, A, f F x(xbr) X A < X br ; , 

^^Um^p 1 (25) 

For the assumed spectral shape, low energy photons (x < Xb r ) are side-scattered at a slightly higher 
rate than are photons near the break energy Xb r (the scattering cross-section is not KN-suppressed) . 
These low energy photons then pair produce off photons of energy x 3> Xb r , and deposit a relatively 
large momentum. It will be convenient to introduce the modified parameter 



T ' = Tc x = aT ^ t ^^L (26) 

F m e cr 



which rescales r c by the ratio of the energy flux at frequency x^ T to the total energy flux. 

The energy of the side-scattered photon cannot, however, be taken to be arbitrarily small. 
When xf <C 1, the radiation shell becomes optically thin to pair creation. This leads to an upper 
bound on the energy of the photo-pair, which we now evaluate. Each side-scattered photon overlaps 
the shell for a time ~ Ai/(1 — cos# s ), but the collision rate is proportional to (1 — cos0 s )a~ rf . The 
optical depth to pair creation therefore depends on 6 S only implicitly through x mi : 



T-yj(x s 6 S ) — of, dx . <J~fj X / x m i n {%s )@s) 



UX A X[ 

m e c< J x A in(x B t0s) 



k 2 (a)F x 



x m'm( x s i ®s) 2 ^ 1' i^) 

J 17l e C z 



The dimensionless function k 2 is plotted against spectral index a in Fig. 1. Defining a characteristic 
optical depth t' c (eq. f26|) through the photon shell at frequency x^ r , and noting that in general 
x A in 3> Xbr, one gets an upper bound on the mean energy x pair = (x A + x B ) ~ (x A ), 



pair 



Xb r 



k 2 — 

Xb 



Jl lV(<«-l) 

(28) 
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This upper bound is approximately independent of the scattering 6 S (even though the most probable 
energy of the side-scattered photon that collides to form a pair of energy x^x does depend on 9 S ). 
The corresponding lower bound on the energy of the scattered photon is 



x? > 



2k 



4kT 2 



J-maxl 



pair 



(29) 



c (l - cos6> s ) 

This energy sits below (above) the break energy x^ when T is less than (greater than) 2T* , where 

1 



4fcV2 



{x^X 



pair \l/2 
max/ 



(30) 



Thus, the supply of side-scattered photons is smaller when V > 2T* , and the acceleration rate is 
reduced accordingly. 

We first calculate the acceleration at low Lorentz factors, r < P, The softness of the side- 



scattered photons allows us to perform the integral (22) over xf in the Thomson regime, dxf /xf 

r Efax ^ xB ■ Setting a lower bound (|29"| ) to x 



dx B /x B . At this stage of the acceleration, (x A ) 
we obtain 



f/3 x max Fx ( x br 



and we have 



>> 



T pair 
x max 



m e 



1-/3 



(r < r*). 



(31) 



(32) 



Xbr / V /' / At 

It is worth recalling the physical origin of the various factors in this expression: fa is the Compton 



afterburn parameter (eq. |17|) which describes the influence of secondary Compton scattering and 
pair creation on the net momentum injected per high energy pair; t' c (eq. ]26|) is a dimensionless 
parameter, proprotional to the energy flux at frequency Xbr an d comparable to the optical depth 
to 7 — 7 collisions through the photon shell; xgf^ (eq. [p8{j) is the maximum energy of the injected 
pair in the lab frame (above which the optical depth to pair creation through the shell is small); 
/i (eq. ||) is the mean inertia per scattering charge; and At is the temporal width of the photon 
shell. 

To approach the high-Lorentz factor regime (where T 3> 2r*) we approximate the scattering 
charges as being cold, ( r y e ) = 1 and \i e = m e . This approximation turns out to be justified, as in 
this regime the time for freshly injected pairs to Compton cool in the bulk frame is short compared 
with the acceleration time (§ |2.2|) . Since the minimum value of x B lies above x^ r , the frequency 
dependence of F x gives an additional factor of 



1 

a 



xP^x br (l-/3)(l-cos fl s ) 
2k(l - (3cos9 s ) 



-i a-l 



(33) 



after performing the integral over xf. Here, the factor in brackets is x^j. divided by the minimum 
value of x B before scattering (eqs. || and f2~9| ). After substituting the relation (1 — cos# s )/(l — 
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Pcos9 s ) = (1 — cos0g)/(l + 0), the acceleration rate can be expressed as 



d 



2fpk 



dt 

which becomes 



a 



2x jjj- /c 



dt 



F x (x hT ) 



,,pair 



(1-/3) / dW 



(1 -cos^ 



1-/3 



2(r*) 2 Ai 

for a = 2. This expression matches onto (32) at T = T*. 



1 + 

(r » r*) 



da 



(34) 



(35) 



Let us compare the rates of acceleration just derived with the corresponding expression 

rR 



dT 

dr 



Y\i-py 



,:2 ■ 



(36) 



for Thomson scattering without pair creation (Noerdlinger 1974; Paper I). When 1 <C T <C T*, the 
acceleration rate (32) is larger by ~ 2fpx^ yi than for Thomson scattering by a spectrum localized 
near x^ T ; and when r 3> T* the acceleration rate (|35"1) is larger by ~ 2fpx^^ x (T/F*)~ 2 . 



Pair creation no longer increases the speed of the bulk flow when two colliding photons of 
comparable energy, x B ~ x A ~ 2k/ xf (1 - cos6» s ) ~ x^^^L have an optical depth to pair creation 
that is close to unity. The resulting maximum Lorentz factor is 



r pair ~ 

max 



X 



pair 



1/2 






k 2 x^ 2 r' c 



l/(a-l) 



(37) 



Continued scattering in the pair- loaded flow increases its speed further: the limiting Lorentz factor 
is 

r™.~-T C) (38) 



from eq. (||). This is comparable to (j3^) for a = 2, but larger for softer high energy spectra. We 
conclude that the main effects of pair creation are i) to increase the rate of acceleration at low 
Lorentz factor; and ii) to reduce the mean mass per scattering charge to fj, ~ m e , thereby allowing 
much higher terminal Lorentz factors. 



2.2. Mean Energy of the Pairs in the Bulk Frame 

It remains to calculate the mean energy {^ e ) of the pairs in the bulk frame, and show that our 
starting assumption that (je) <C is self-consistent. During the first stages of acceleration an 
electron or positron will be injected with a characteristic energy <-y™ ax ~ 2 x n5tx> an d w ih co °l off 
the radial photon flux when x^ r <C ^max- At the same time, the accelerated medium is compressed 
(Paper I), and the pairs are adiabatically heated at a rate {^ e )~ 1 (d{-f e ) /dt) = \n'~ x (dn' / dt) = 
^T~ 1 (dT/dt). The equilibrium energy (j e ) that results from a balance between compressional 
heating and Compton cooling lies far below the injection energy, and so we will assume that the 
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heating process has isotropized the momenta of the pairs in the bulk frame. Prom eq. (38) of Paper 
I, 

1 ft V /2 

<7e> = -z hp 1 , (39) 

Z \ ^accel / 

where the reference (lab-frame) cooling time is ijL j = rrieC 2 / utYFy ~ (r/2r^)At, and Fr — F/2Y 2 
is the energy flux in the boosted frame. The acceleration rate is decreased by a factor ~ (7e) 1 
from the values (|32|) and (^) calculated above for a cold plasma: 



1 dY 1 



Ydt ( 7e )t° 



accel 



For a spectrum a = 2, 



atl«r«r, and 



E accel ^ X br Vi 



Ci _*//?* (Y 



f „2 \ r * 

''accel x br vi 



(40) 



(41) 



(42) 



at T > r*. The equilibrium Lorentz factor works out to (je) ~ 1.4T 1 (x^Hr/WO) 1 ^ 2 for a pair- 
loaded plasma moving slowly (Y < r*). 



2.3. Relative Inertia in Pairs and Baryons 

The distribution of inertia between (cold) pairs and baryons within the accelerated material 
has important implications for the radiative efficiency of the forward shock. The number of pairs 
multiplies greatly as each hard (of bulk frame energy j e ~ x^^/AY) cools by side-scattering 
the photon beam. During the initial stages of the acceleration, the side-scattered photons are well 
above the pair creation threshold; little additional momentum is absorbed from the photon beam 
when they collide with soft photons. The momentum imparted to the medium by the cooling of 
the hard is cancelled. In this regime, / — 1. 

Balancing the radial momentum density in the lab frame, under the assumption that each 
injected hard pair leads to iV mu it 2> 1 softer pairs with an isotropic distribution in the bulk frame, 
we obtain 

^max -TF^— meC = r ( 2 ^n e + + m p n p ) c. (43) 

Here, \x e = |(7 e )m e when the pairs are at least mildly relativistic and the multiplicity iV mu it is 
given by a;^w/4r(7 e ) (eq. [fl3"l] ). This equation inverts to 

= 2. (44) 
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The medium becomes very heavily loaded by pairs as it accelerates to bulk relativistic speed (where 
{j e ) ~ 1). The fraction of the inertia carried by the hadrons becomes 



& = . m 'l' A («) 



in this regime. 



Compton scattering of photons of lab-frame energy ~ x^ r by the hard is initially in the 
Klein-Nishina regime, where the loses a significant fraction of its energy in one scattering. 
Cooling continues in the KN regime as long as 7g(xb r /2r) > 7 e , that is, as long as T < (2k) l l 2 T* . 
As the medium attains a Lorentz factor not far below the maximum value (p8|), the hard are 
injected with a small enough bulk frame energy that most of the energy lost to Compton cooling 
is not reabsorbed by photon collisions. This occurs when (ar te /2r)(xg£/2T) < 1, that is, when 
(r/r max ) 4 > \ k 3 kl /{a ~ 1] (r^/xbr)^- 40 )/^- 1 ). Above this Lorentz factor, / ~ f . 

Implicit in this calculation of the maximum Lorentz factor is the assumption that the photon 
flux is high enough to supply the needed momentum. This constraint translates into an upper 
bound on the density of the medium surrounding the 7-ray source. The strongest constraint comes 
from the requirement that the phase of pair-dominated acceleration continues up to a Lorentz factor 
TPfax (eq. pT^I ). This sets a lower bound to the flux of photons of energy ~ a^ax- 



3. Shock Deceleration, Pulse Broadening, and Afterglow 

We now turn to consider the effects of radiative acceleration on the dynamics of a decelerating, 
relativistic fireball. The emission process operating in a GRB probably covers a range of radius. In 
some models, the emission is due to optically thin synchrotron and/or inverse-Compton processes 
(Meszaros, Rees, & Papathanassiou 1994; Sari, Piran, & Narayan 1998); whereas in others the 
emission is by inverse-Compton scattering at moderate Compton optical depth (Thompson 1994, 
1997; Crider et al. 1997). To focus the present discussion, we assume that the dissipation extends 
outward from (at least moderately) large optical depth. The optical depth r c (eq. |§) through the 
photon shell is 

Tc = T?(i-Y\ (46) 



outside the scattering photosphere. Here, R T is the radius at which the ejecta become optically 
thin to scattering and the 7-ray photons begin to stream ahead of the matter. As long as pair 
creation and annihilation remain in equilibrium at energy x ~ 1 in the bulk frame (x rs -' r e j in the 
lab frame), the position of the Thomson photosphere is determined implicitly by^ 



Fjr ei x h y-^ = ^(^) ~1. (47) 



Here a is the high-energy spectral index; eq. (b 
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(Recall that t' c is r c rescaled by the factor x\yrF Xh /F.) In this way, the initial value of t c can be 
related to the (initial) bulk Lorentz factor of the ejecta, 

^"(Tj) - 1 *^. (48) 

This equation provides an upper bound on if the 7-ray emission region is optically thin to 
scattering. 

The width of the prompt 7-ray pulse is bounded below by 



A4 >- ws?- (49> 



More generally, we can define a minimum timescale for variability driven by internal dissipation in 
the outflow, 



^r{R) = 2^2;. (50) 

This timescale can be related to the 7-ray fluence dE/dVt (at energy Xb r and per unit solid angle) 
by noting that 



dE/dfl 



< = vt^jj^- (51) 
At the scattering photosphere, one finds 



At var (i? r ) = 0.5 . , 



10 53 erg 



1/2 / F \ -(3+a)/2 



300 



(52) 



upon making use of eq. fl4q)- Should the overall duration At of the burst be determined by the 
activity of the central engine, one does not expect the afterglow emission to connect up smoothly 
with the main burst. However, the internal dissipation could be triggered by an interaction with 
the external medium: for example, when a blob of magnetized ejecta is decelerated and the causal 
propagation distance within the blob grows large enough to allow reconnection (Thompson 1994; 
Begelman 1998). In such a hybrid model, the observation of distinct, overlapping sub-pulses requires 
that the ejecta extend over a cone of angular width 9 3> (r^) -1 , and that the surface of the ejecta is 
wrinkled on angular scales larger than ~ (r[!j) -1 . To simply the discussion that follows, we assume 

,f width At ~ R_/9.fr - 



that the burst is composed of a single pulse of width At ~ R T /2(T^-) 2 c 



3.1. Deceleration of the Ejecta 

Pre-acceleration of the ambient medium by the prompt 7-ray pulse will slow the deceleration 
of the burst ejecta, and therefore modify the synchroton emission from the forward shock. The 
maximum Lorentz factor (eq. p8| ) to which the ambient medium may be accelerated decreases 
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with distance from the matter photosphere. Let us assume, for the time being, that the ambient 
material is light enough to allow acceleration to r am b = r max . Then 

r amb («) = r max (i? r ) (4-) , (53) 



. Rt 

since r max oc r c oc R~ 2 . Thus, the kinetic energy per scattering charge behind the forward shock is 

-LLC, (54) 



Tamb(^) 

in the lab frame. Here, fj, = m e /(l — fh) (eq. ||) is the mean mass per scattering charge. 

The energy E of the shock decreases with radius, because only a fraction fraction fh of the 
post-shock particle energy is carried by baryons that do not cool. The calculation in §2.3 indicates 
that the inertia of the ambient medium becomes dominated by leptons (eq. ^) . Energy is deposited 
in shocked leptons at a rate 

S- (1 - A) ra"- (fi)fl2c3 (55> 

We assume that these particless cool instantly. When fh <C 1, the remaining energy E of the 
fireball is dominated by the kinetic energy of the ejected matter (mass Mo), 

^-r-^V (56) 



One can calculate the time-scaling of the shock energy by balancing the energy radiated ( |55j ) against 
the energy retained (56). Assuming the density profile around the center of the explosion to be a 
powerlaw 

n p0 (R) = ^ oc R- 5 , (57) 

one finds 

E{R) oc R-% (58) 

with 

e = 5 - 5. (59) 

This compares with e ~ 3 — d in the case where the external medium is static or moves subrela- 
tivistically (e.g. Katz & Piran 1997). 



More generally we will assume a power-law scaling ( pq ) of the shock energy. This assumption 
is motivated by the weak dependence of fh (eq. |^5|]) on the gamma-ray continuum in the region 
where acceleration to bulk relativistic speed is possible. At large radius fh — ► 1 and e — > 0, because 
pairs with the miminum (shock-frame) energy ~ F e ^m e c 2 no longer cool effectively. This allows a 
first evaluation of the influence of radiative acceleration and pair loading on the initial stages of 
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the deceleration of a relativistic shock. A more detailed numerical treatment of shock deceleration 
including these effects is deferred to a later paper. 

Balancing the shock energy with the energy deposited in shocked particles, one obtains 

r ej (*) = r,(^) . (60) 

The relation between shock radius and observed (post-burst) time t becomes 

t oc R 6 - s+t . (61) 

Pair creation feeds back strongly on the radiative losses from a relativistic outflow, and regions 
with extended non-thermal spectra are expected to be much brighter than regions with quasi- 
thermal (e.g. Wien) spectra (Thompson 1997). It also introduces new characteristic lengthscales 
and time-scales into the afterglow process. One of these is the radius -R(r am b ~ 1) outside of which 
the ambient medium is accelerated to sub-relativistic speeds. If the medium were to remain pair 
loaded out to this radius, one would calculate -R(r am b ~ 1) by setting [i ~ m e in eq. ([!]); this 
would yield i?(r am b ~ 1)/R T ~ {\ T c) 1 ^ 2 - However, because several efoldings of the pair density 
are required to push the ambient medium (assumed initially free of pairs) to bulk relativistic 
speed, bulk relativistic motion effectively stops somewhat inside that radius. A more accurate 
estimate of -R(r am b ~ 1) is obtained by noting that when r c is not much larger than unity, the 
pair density exponentiates according to n e +/n p ~ exp(n 7 o"TcAi) = exp(r c /xb r )- Requiring that 
2m e n e + ~ mpn p yields r c /xbr ~ ln(m p /2m e ) ~ 7. This value of r c is large enough that the photon 
shell is optically thick to photon collisions, r 77 ~ 0.25T c /xb r ~ 2. Pairs continue to dominate 
baryons in number out to a radius where r c /xb r ~ 1; but are not produced in sufficient numbers 
to force the mean inertia per scattering charge down to /i ~ m e . The corresponding delay from 
the beginning of the burst is given by the time it takes the forward shock to reach that radius. 
In the case of a uniform external medium, T e j oc R~ 5 ^ 2 and the time delay scales with radius as 
R/2Tl-cocR 6 . This yields 

A t(r b ~i) /^y 

At \ 10 / I 300 / K J 



3.2. Synchrotron Emission 

This pre-acceleration of the ambient medium will modify the time-scaling of the synchrotron 
emission from the forward shock propagating ahead of the burst ejecta. The pair-loading ensures 
that the Lorentz factor of the shock-accelerated particles (in the shock frame) extends downward 
to the minimum kinematically allowed value, j e ~ r e j (Thompson 1997). We consider here only 
the simplest emission model involving synchrotron radiation by a power-law distribution of pairs, 
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dn e /dj e oc 7 e for rest frame energies 

7e > (63) 

1 amb 

Closely related synchrotron models, in which the inertia of the non-thermal particles is dominated 
by baryons, have been constructed by Meszaros k, Rees (1997), Wijers, Rees & Meszaros (1997) 
and Waxman (1997). The energy carried by the non-thermal pairs is taken to be a constant fraction 
e cr of the total post-shock particle energy in the shock frame, and the magnetic field is assumed to 
carry a constant fraction eb of the post-shock pressure. As we discuss below, pre-acceleration of the 
ambient medium by a collimated 7-ray beam will ensure that this second assumption is satisfied, 
once the ejecta are moving sufficiently slowly. 

The medium just ahead of the shock is compressed at the same time as it is accelerated, so 
that the lab-frame density of baryonic particles has increased to n p ~ 2r^ mb n p o over the ambient 
value (eq. [17] of paper I). The bulk frame energy density behind the shock is, then 

(fr^-) 2n e+/ uc 2 ~ — ^-n p0 m p c 2 . (64) 

V J- amb / Jh 

The equilibrium post-shock magnetic pressure therefore scales as 

B 2 

g£ oc e B T%n p() oc iT 5 " 6 , (65) 

independent of S in the shock rest frame. 

We are now in a position to work out the time-scaling of the synchrotron emission at fixed 
frequency v. Pairs which radiate at this frequency (Lorentz-boosted into the observer's frame) have 
an energy 



7e 



/2^m eC \ 1/2 o(i?5/2 _ (5 _ 2e)/4 ^ (66) 



When synchrotron cooling is rapid, the synchrotron emissivity depends on the incoming flux of 
kinetic energy, but not explicitly on post-shock magnetic field. Assuming an powerlaw spectrum 
dN e /d"/ e oc 7~ p of injected particles, the energy released in photons of frequency v is, per unit 
time 



R 2 on ~ CJ 



2 



\ 2 ~ P 3 
Jl 1 or2^ <-D\ m v c 



2r e > p0 ( J R)- 



(67) 



1 \ r /r ; e J v y ' f 

\ 1 ej / 1 amb / J h 

The quantity in brackets denotes the rate at which particle energy accumulates (in the bulk frame) 
per unit area of the shock; the prefactor transforms to the lab frame. The power law index 

is expressed in terms of the synchrotron index a sync = \P via 
6a sync + 2 - i(5(3a sync - 1) + 2ea sync 

P = — R F" — (r a mb = r max ; fast cooling). (69) 

6 — + e 
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When the shock loses energy slowly (e ~ 0), one has (3 = a sync + | for a uniform medium (5 = 0) 
and (3 = |a sy nc + | for a steady wind (6 = 2). 

Now let us compare this result for the time-scaling of the synchrotron emission, with the 
corresponding expression for a shock propagating into a static external medium. In that more 
familiar case, the shock Lorentz factor, observer's time, and rest-frame magnetic field scale as 

r ej OC R^-^' 2 , t OC i?V(4-<5+e) ) Br K #-3/2-e/2 ; &nd ^ K R 3/2-(S-2e)/4 _ For thoge sca l ingS) foe 

time index is 

a 6a sync - 2 - ^(3a syn c - 1) + 2ea syn c , , . . 

P = 4_^ + e ( ramb = 1 ' fast coolm g)- ( 70 ) 

This agrees with the scaling previously obtained whose energy decreases by Sari et al. (1998) and 
Meszaros et al. (1998) if we substitute expression (|5^ ) for e (in the case of a strongly radiative 
shock); or e = (in the case of an almost adiabatic shock). An adiabatic shock has (3 = §a syn c — \ 
for any value of 5. 

The regime of slow synchrotron cooling (at the observed frequency v) is treated analogously. 
We can now assume e = 0, because the lowest energy must not be able to cool in this regime. 
One has 

dL v ( d 2 N e \ 2 2 2 
where the total number of radiating charges is 

a ^ fl Vo(A) ^ (rfr-) l " ■ (72) 

d-y e dil 3-d fame \T e] /T amh J 

The relation between synchrotron index and particle index softens to a S ync = (P - l)/2, and the 
time index becomes identical to (|69|): 

6a sync + 2 - i(5(3a sync - 1) 
/3 = y - 2 (r amb = T max ; slow cooling) (73) 

D — 

for e = 0. Notice that the functional dependence of the time-index (3 on the spectral index a sync 
remains the same for both fast and slow synchrotron cooling, when the ambient medium has been 
pre-accelerated by the prompt 7-ray pulse and the lowest energy e ± do not cool radiatively (e ~ 0) . 
If the injected spectrum of non-thermal particles is constant - let us take P = 2 - then across a 
'cooling break' the spectrum hardens from ct sync = 1 to Qsync — 5 an d the time index hardens from 
/3 = (8 - <5)/(6 - 5) to (5 - \8)/{<o - 8). 

By contrast, there is a stronger break in the time scaling between regimes of fast and slow 
cooling if the external medium is static, 

P = 6 " SynC " f^ SyDC " 1} (r amb = l; slow cooling). (74) 

(Meszaros et al. 1998). This simplifies to the familiar result /3 = |a sync in a uniform external 
medium. For P = 2, the time index hardens from (3 = 1 to (3 = (3 — \S)/{A — S) across the 'cooling 
break'. 
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3.3. Anisotropic Radiation Pressure and Shearing of External Magnetic Fields 

The prompt 7-ray pulse induces strong shearing motions in the ambient medium that can 
amplify a seed magnetic field before the forward shock hits. An important question, that has yet to 
be resolved in afterglow models, regards the effectiveness with which a weak ambient magnetic field 
(e.g. B ex ~ 3fi G for a fireball interacting with the ISM) will be amplified above the flux density 
~ T^-B ex expected from laminar compression behind a relativistic shock. In the non-relativistic limit 
(where pair creation can be neglected), the limiting speed of the ambient medium is proportional 
to the 7-ray flux. One therefore expects the flow outside the forward shock to be strongly sheared 
on lengthscales in between (rgj) _1 i? and #beam-^! given that the 7-ray emission is beamed within 
an angle #beam but also varies on the smaller angular scale due to causal fluctuations in the rate of 
internal dissipation. This sheared flow will be Kelvin-Helmholtz unstable and, if there is time before 
the shock hits, the resulting turbulence will strongly tangle a seed magnetic field. Since the available 
time is ~ R/2T^c, one expects effective tangling only after the ejecta have decelerated below 
T c j ~ (r^j) 1 / 2 . Thereafter, near equipartition between the magnetic pressure and the turbulent 
pressure in the sheared flow appears plausible: B ~ (AnnpompC 2 ) 1 ^ 2 ~ 0.1(n p o/l cm -3 ) 1 / 2 G. If 
no other instability (such as a two-stream instability; Medvedev & Loeb 1999) raises the field close 
to equipartition, one expects the burst and/or afterglow lightcurves to contain features correlated 
with the onset of effective magnetic shearing. 



3.4. Upper Bound on Pre-Burst Mass Loss: Pulse Broadening 

A strong constraint on the density of the ambient medium comes from the requirement that the 
accelerated matter (which becomes loaded with pairs) remain optically thin to Compton scattering, 
so that the 7-ray pulse is not appreciably broadened. This becomes a serious constraint even if the 
scattering optical depth through the ambient medium is small (cf. §3.3 of Paper I). If all photons 
of frequency > x are converted to pairs, then the scattering depth through the resulting shell of 
pairs is 

tL ( x V~ a tL° ( R \- 2 ( x ^ l ~ a 



TT ^^L _ = jl. _ . (75) 

Xbr \XhrJ X hr \R T J \X hr J 

In order to prevent photons of this energy from being completely consumed, the medium must be 
accelerated sufficiently that the photons are no longer above the threshold energy for pair-creation 
in the bulk frame.^] To maintain tt < 1, we require that 

r amb > r T (R) = \t' c {R) x«; 2 ] 1/{a ~ 1] . (76) 



5 The material undergoing radiative acceleration remains in causal contact with itself, and so the scattering depth 
through it is not affected by the suppression of the scattering rate (by a factor ~ 1 — f3) in the lab frame. 
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The critical Lorentz factor T T can be expressed more directly in terms of the initial bulk Lorentz 
factor of the ejecta, after making use of eq. (48): 

r (r) ( r \ -Vfa- 1 ) 

tWJ ' . (77) 



Note that T T coincides with at the scattering photosphere. If T am b < IV, then the 7-ray pulse 
is spread out to a width 

At spre ad _ R /r amb \- 2 / fl\(a+3)/(a-l) 

At 2r2 mb cAt "IrJ liij • U8j 

One finds that At sprea d > At at R > R T . 

The kinetic energy of the accelerated medium moving at Lorentz factor IV (within a shell of 
thickness equal to the acceleration length fT^cAi; eq. P]) can be no larger than the energy of the 
7-ray pulse above a frequency x ~ IV, 

f- me] 



r T n p0 m p c 2 f 1 + 2 ^i^ j -r 2 cAt < f^^) (IlV ° c At. (79) 



This yields the upper bound 

o (5-2o)/(a-l) / . / n + m \ _1 /R\6/(a-l) 

«<%iF^r(^)( 1 + 2 ^) K (D ' (80> 

The density in the above equation is higher than that encountered in the diffuse ISM, and therefore 
does not provide any constraint on GRB models involving the merging of compact objects. We can 
use this expression, however, to set limits on massive stars as likely GRB progenitors, as in this 
case the explosion occurs in a dense pre-burst stellar wind. Assuming a steady mass loss from such 
a massive progenitor at a rate M and constant speed V, the ambient density is 



Substituting this expression into eq. (pQ) , and making use of eqs. ( |4q) and (|49j) to express R T in 
terms of At, t' c ° and xt, r , we have 



/ \ l/(a-l) m rVAt / R \ (2<*+4)/ (a— 1) 

M<24vrA(x 2 r 3 r^) 1/(Q (if) 



zj 



(82) 



Setting T am b = T r and eliminating , one deduces that the amount of pre-burst mass loss needed 
to induce a certain spreading of the pulse is 



M = / At spread \ 
M A t \ At ) 



(2a+4)/(3+a) 



(83) 
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where 



M At = 247rf h x hr T° ei 



in, r\ Ai^ 



1.5 x 10- 5 x hl f h ) ( —) ( - r ) M yr- 1 . (84) 

J 1 300 / VlOs/ V 1000 km s- 1 ' 



Notice that both the width At sprea d of the spread pulse (eq. |49|| ) and Ma* grow with radius R. 
Mass loss rates close to the critical value Ma* only broaden the 7-ray pulse near the scattering 
photosphere, and do not broaden it by much. Higher mass loss rates cause further broadening at 
larger radii. 

The presence of a dense medium surrounding the GRB source also limits the bulk Lorentz 
factor Tgj of the ejecta at the Thomson photosphere. Hence the width At of the emitted 7-ray 
pulse (before scattering in the ambient medium) cannot be made arbitrarily small. An upper 
bound to Tgj is obtained by balancing the energy of the ejecta with the kinetic energy of the 
baryons accumulating behind the shock: 

' ^>^Po(^c 2 . (85) 



e rad 89, 3-5' 

Here e ra d is the fraction of the kinetic energy of the ejecta converted to 7-rays; this expression is 
meant to be evaluated at the scattering photosphere of the outflow, and so is missing the factor of 
Tamb present in eq. |^6| . Combining this equation with (|48|) and ( 51 ) , one finds 

( j^v_ ) v<mo /aw \ 



or 



T°: < 36 e 2/3 



2/3 / V \ 2 / 3 



M \ ( V 



ej - ^ c rad i n _ 5 A/ r „_l 



10- 5 M Q yr- 1 I V 1000 km s 



4ir(dE/dn) 



10 53 erg 



1/3 



(87) 



for a hard spectrum a = 2. The width of the broadened pulse becomes 

^ ~ 2(r°.) 2 c ~ 2c { AttV ) \m) \o T ^ ) 

Evaluating this expression for a = 2, we deduce 



At > 170 e 



5/3 I M V /3 ( V V 5/3 \M9E/dQ) y 1/3 



rad I 10 _ 5 Mq yr _i I yiQOO km S" 1 



10 53 erg 



s. 



We can combine this last result with eq. fl83|) to obtain a very similar bound on the width of the 
scattered 7-ray pulse: 

^spread = ( 1 \ (a + 3)(a-2)/(2a+4)(5- a ) . , 

At v3Ax br e rad ; ^ ■ K ' 
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Note that the width of the scattered pulse is independent of T e j for a hard spectrum a = 2, 
At sprea d/Ai = 2.2 (A^bAad/O-l) -5 / 8 - It is much larger than At for softer spectra, e.g., Ai spreat i/Ai 
= 2.1 T^ 10 (A^br^rad/O-i) -3 / 5 for a = 3. Because the softer spectrum creates fewer pairs, a larger 
external mass density is required to induce runaway pair creation for a fixed shock energy. 

This calculation sets significant constraints on GRB models involving the delayed explosions of 
massive Wolf-Rayet stars (Woosley 1993; MacFadyen & Woosley 1998), whose winds are character- 
ized by mass loss rates M 1(T 5 - 1CT 4 M & yr~ x and velocities of 1000-2500 km s' 1 (Willis 1991). 
Hard bursts composed of peaks narrower than ~ 10 sec must either form in a more rarefied envi- 
ronment, or have a low radiative efficiency e ra d- The possibility remains that the overall duration 
of some 7-ray bursts is determined by the processes of pair-loading and radiative acceleration in 
the ambient medium. We note finally that the constraint (|89|), ( |90|) is weakened if the 7-ray source 
function is not a rigid power-law, but has a thermal cutoff below a bulk frame energy ~ m e c 2 . 

4. Discussion 

We have calculated the deposition of radiative momentum into an (initially) static medium 
by the intense 7-ray pulse emanating from a relativistic fireball (eqs. |32| and |35|1 ). This process 
loads the accelerated medium with pairs, which come to dominate its inertia (eq. ^4|); and greatly 
increases the radiative efficiency of the forward shock wave that bounds the fireball material. Given 
this high radiative efficiency, the ambient medium will attain a bulk Lorentz factor close to the 
Lorentz factor of the fireball material near the 7-ray photosphere, decreasing as R~ 2 at larger radius 
(eq. |p3[ ). As a result, the rate of deceleration of the ejecta is slowed (eq. [p0[]), and the relation 
between the time-index of the synchrotron emission and the synchrotron spectral index (eqs. 
and f73|) is modified from the relation previously obtained for a static ambient medium. 

These physical processes have important implications for the radiative physics of the forward 
shock, aside from its overall radiative efficiency. Anisotropy in the 7-ray flux (on an angular scale 
r _1 or larger) will induce shearing motions in the external medium that can strongly amplify 
any seed magnetic field before the shock hits it (§ |3.3j ). An even more important effect is that the 
minimum Lorentz factor of the shocked pairs is pushed down to 7 e ~ 1 in the bulk frame, so that the 
corresponding synchrotron break frequency lies at optical-UV near the e^ scattering photosphere. 
This may point to a hybrid radiative mechanism (Thompson 1997) in which hard X-ray photons 
advected out with the fireball are the principal seeds for Comptonization by non-thermal pairs and 
bulk fluid motions. 

Support for this work was provided by the Sloan foundation (C. T.), and by the NSF through 
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Fig. 1. — Left: The function k(a) (eq. 24]), plotted against spectral index a. Right: The function 
k 2 (a) (eq. 



